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ABSTRACT 

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in 
a forming cluster to those in the field at z = 1.6 with Hubble Space Telescope near-infrared imaging in the 
CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, r e ff, and Sersic 
index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (« > 2) reside on the 
red sequence of the color-magnitude diagram, and most disk- dominated galaxies (« < 2) have colors expected 
for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U —B colors 
and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no 
significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster 
have larger median effective radii and smaller Sersic indices compared to the field with a significance of 2er. 
These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < R pm j < 1 .5 Mpc, which 
have low Sersic indices and possibly larger effective radii, more consistent with star-forming galaxies at this 
epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under 
the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are 
consistent with models where gas accretion onto these galaxies is suppressed from processes associated with 
the cluster environment. 
Subject headings: galaxies: high-redshift 
structure 
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1. INTRODUCTION 

Galaxies in the low redshift Universe show a very strong 
dependence between environment, morphology, color, and 
stellar mass (e.g., Hogg et al. 2004; Kauffmann et al. 2004; 
Skibba et al. 2009). Galaxies in regions of higher local den- 
sity, such as galaxy groups and clusters, have lower levels of 
star-formation, higher fractions of red, passive galaxies, and 
more early-type morphologies compared to galaxies in lower 
local number density regions (e.g., Dressier 1980; Postman 
& Geller 1984). Poggianti et al. (2008, and others) find that 
this relation extends to intermediate redshifts (z ~ 0.4-0.8) 
with the highest density regions being exclusively populated 
by elliptical morphological types. 

It is difficult, however, to separate the galaxy properties 
which depend on high local density (regions of dense "en- 
vironment") from those which depend on galaxy mass. One 
challenge is that galaxies in regions of high local density 
have higher average halo and stellar masses, and care must 



be taken to study environmental effects on samples at fixed 
stellar mass. For example, Kauffmann et al. (2004) showed 
that the red, passive galaxies at z ~ 0.1 that dominate in high 
local density regions have highly concentrated morphologies, 
but at fixed mass their sizes and concentrations are almost in- 
dependent of environment. Guo et al. (2009) find that the 
concentrations of the central galaxies of dark matter halos at 
< z < 0.2 depends strongly on their stellar mass, but only 
weakly on their halo mass, implying a weak dependence on 
environment. Similarly, van der Wei (2008) show that galax- 
ies at z ~ 0.02 with the highest masses show no dependence 
between morphological concentration and environment (lo- 
cal density of galaxies). It is only for lower mass galaxies 
(M~3 x 10 10 -5 x 10 10 M Q ) that there is a weak trend toward 
higher galaxy concentration in higher local density regions. 
Therefore, the environment has a subtle effect on galaxy mor- 
phologies at low redshifts. 

Griitzbauch et al. (2011a) show that these trends persist to 
higher redshifts (z ~ 0.4- 1 .0) and that the intrinsic properties 
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of a galaxy are more dependent on stellar mass than local den- 
sity at earlier epochs. In a similar redshift range some stud- 
ies attempted to disentangle the color-mass-density relation 
and its evolution by examining the colors or morphologies of 
galaxies as a function of local density in narrow bins of red- 
shift and stellar mass (Tasca et al. 2009; Cucciati et al. 2010). 
These studies found that the correlation between color, mass, 
and density is considerably weaker at higher redshifts with the 
exception of galaxies in the low mass regime. The implication 
is that these lower mass galaxies were formed more recently 
in a time when evolved large-scale structures were already in 
place. 

In contrast, the environment does appear to have a strong 
effect on the relative star-formation rates of galaxies. Hogg 
et al. (2003) showed that at z ~ 0.2 the density of galaxies cor- 
relates with color: redder galaxies lie in regions of higher lo- 
cal density (see also, Blanton et al. 2005; Baldry et al. 2006). 
Similarly, Kauffmann et al. (2004) found that higher density 
regions anticorrelate most strongly with galaxy specific SFR 
(sSFR, the SFR per unit stellar mass): galaxies in high local 
number density regions have lower average sSFRs. Patel et al. 
(2009) show that this trend also persists to higher redshifts 
(z ~ 0.7). Interestingly, most of this effect on the sSFRs does 
not arise from an environmental dependence on the properties 
of star-forming galaxies, which appear unchanged between 
regions of high and low local number density (Hansen et al. 
2009; Peng et al. 2010; Wetzel et al. 2011). Instead, it is the 
fraction of quenched galaxies that increases with increasing 
galaxy density that drives this observation (Peng et al. 2010). 

Several authors have argued that this quenching fraction is 
different for galaxies that are "centrals" or "satellites" in their 
halos, where quenching in centrals depends primarily on the 
galaxy's halo mass while the quenched fraction of galaxies 
that are "satellites" depends on their distance from the halo 
center (e.g., Wetzel et al. 2011; Wilman et al. 2011; Woo 
et al. 2012). Hansen et al. (2009) find that the fraction of red 
(quenched) and blue (star-forming) galaxies changes dramat- 
ically with distance from the center of clusters. Wetzel et al. 
(201 1) find that the fraction of quenched galaxies persists out 
to ps 2R WU , where R wn is the virial radius of the halo (e.g., 
Hansen et al. 2009). Selecting galaxies in bins of morpholog- 
ical type, Guo et al. (2009) find no evidence that the struc- 
tures of early-type galaxies differ between central and satel- 
lite galaxies matched in both optical color and stellar mass, 
although Weinmann et al. (2009) find that late-type satellite 
galaxies show smaller radii and larger concentrations than 
late-type centrals. Weinmann et al. (2009) also find that late- 
type satellites have redder colors than late-type central galax- 
ies. This suggests that the stellar mass of the galaxy is the 
more fundamental property determining the structure of early- 
type galaxies, and that observed differences in satellite and 
central galaxies may be due to the quenching of galaxies after 
they become satellites. Woo et al. (2012) argue these observa- 
tions are consistent with theoretical expectations that quench- 
ing mechanisms affect the satellites as they move into the en- 
vironment of the larger mass halo (e.g., Bahe et al. 2012). 

It is unclear how these trends between color, SFR, morphol- 
ogy, mass and environment evolve with redshift. Both Cooper 
et al. (2007) and Elbaz et al. (2007) observed that the low 
redshift SFR-density relation evolves and possibly reverses 
around z ~ 1 . It has also been observed that at z ~ 1-1.6 
the surface density of IR-luminous, star-forming galaxies in- 
creases in regions of high local density (e.g., Tran et al. 2010; 
Kocevski et al. 2011). However, the study of Grutzbauch 



et al. (201 lb) found that in the redshift range 1 < z < 3 there 
is a strong dependence between SFR (and sSFR) and stellar 
mass, and very weak (if any) dependence between SFR and 
environment (see also, Bauer et al. 2011a). It is only in the 
highest local number density regions where star formation is 
significantly suppressed, and this is only apparent to z ~ 2. 
However, studies show that even at these redshifts, regions of 
high local density contain a dominant population of mostly 
quiescent galaxies, and the overall trend in color and specific 
SFR mostly unchanged from z ~ (Cooper et al. 2008; Chuter 
et al. 201 1; Quadri et al. 2012). 

Other studies have looked at the relation between galaxy 
size (morphology) and environment at z > 1 (Cooper et al. 
2011; Zirm et al. 2011; Papovich et al. 2012). These studies 
have found that at fixed mass, spheroidal (bulge-dominated, 
with Sersic index n > 2) galaxies in regions of higher den- 
sity have larger sizes compared to spheroidal galaxies in low- 
density regions. This seems to imply that there are processes 
associated with the higher density regions that enhance or ac- 
celerate galaxy growth. One obvious candidate is an appar- 
ent higher merger rate for galaxies in higher density environ- 
ments. Indeed, Lotz et al. (2013) observe an enhancement 
in the rate of dissipationless (gas poor, or "dry") mergers in 
galaxies in a forming cluster at z = 1 .62 compared to similarly 
selected galaxies in the field (similar to the findings of Tran 
et al. 2005, at z = 0.83). If these "dry" mergers are in fact the 
dominant driver of the size evolution of early type galaxies, 
they may also affect the galaxy morphologies (e.g., Navarro 
1990), and we may expect to find variation in galaxy structure 
at this epoch. 

One reason it is difficult to disentangle galaxy evolution ef- 
fects arising from "environment" from those which depend on 
"mass", is that different methods are used to define "environ- 
ment" and these are sensitive to the effects of environments on 
different scales. Many methods have appeared in the litera- 
ture such as fixed aperture and nearest neighbor distances, but 
it has been shown that there is considerable scatter between 
them. Muldrew et al. (2012) compare and contrast differ- 
ent environment estimators using mock galaxy samples taken 
from the Millennium simulation. They found that different 
methods probe different aspects of the dark matter halos that 
contain groups and clusters. For example, nearest neighbor 
statistics better study the internal halo properties and effects 
of smaller structures such as group sized halos. The fixed 
aperture method, on the other hand, is found to be a better 
probe of the halo as a whole. However, all of these environ- 
ment statistics are sensitive to errors in galaxy distances (in 
particular redshift errors), which can be the limiting factors in 
this analysis. 

Here, we extend these studies to understand how galaxy 
color, stellar mass, and morphology depend on environment 
at z = 1.6 using a portion of the data from the UKIRT Infrared 
Deep Sky Survey-Ultra-Deep Survey (UKIDSS-UDS) field 
covered with the Hubble Space Telescope (HST) as part of the 
Cosmic Assembly Near-infrared Deep Extragalactic Legacy 
Survey (CANDELS, Grogin et al. 2011; Koekemoer et al. 
2011). These data cover one of the largest contiguous fields 
with WFC3 near-IR imaging. At z = 1 .6 this field covers both 
a forming cluster at z = 1 .62 and galaxies in the field. There- 
fore, we are able to study galaxies within a single redshift slice 
in this large field. This allows us to test effects associated with 
environment over a very large range in local galaxy density, 
from the rich environment of a forming cluster to the low den- 
sity field using a highly homogeneous dataset. Furthermore, 



The Environmental Dependence of the Color-Mass-Morphology Relation at z = 16 



3 



because we study galaxies from a single dataset with uniform 
analysis, any bias from systematics affects all galaxies in a 
similar way. This means we can make robust (relative) com- 
parisons between galaxies in different environments. 

Throughout this paper we report magnitudes measured rela- 
tive to the AB system (Oke & Gunn 1983) and we assume cos- 
mological parameters Cl m = 0.3, fl\ = 0.7, and Ho = 70 km s -1 
Mpc -1 , which give an angular diameter distance of 0.508 Mpc 
arcmin" 1 (8.47 kpc arcsec" 1 ) for z = 1.6. All projected dis- 
tances are given in physical (proper) distances at z = 1.6. 

2. DATA 

For this paper, we use available datasets that cover the 
CANDELS UDS field. The available data in this field in- 
cludes //S7YWFC3 imaging in the F125W and F160W filters 
from CANDELS over 9 'A x 22 ?0 (Grogin et al. 201 1; Koeke- 
moer et al. 201 1), corresponding to a projected, physical area 
4.8 Mpc x 1 1 .2 Mpc at z = 1 .6. 

Here we used a photometric catalog containing photometry 
from the UKIDSS-UDS from Williams et al. (2009) cover- 
ing BVRIz from Subaru/SuprimeCam, JK from UKIDSS, and 
deep Spitzer/IRAC imaging from the Spitzer UDS covering 
3.6-8.0 /im. Our catalog is identical to that used in Papovich 
et al. (2010). In summary, these catalogs reach 5<r limiting 
magnitudes of B < 27.7, R<21.l,i< 26.8, z < 25.5, J < 23.9 
and K < 23.6 mag in apertures of 1 "75 diameter. The catalogs 
are selected in the A"-band with colors measured in fixed aper- 
tures and they include unique aperture corrections for each 
source defined as the difference between the fixed aperture 
and total magnitude derived from the K-band for that source 
(see Williams et al. 2009). We used IRAC catalogs with pho- 
tometry measured in 4" diameter apertures, with aperture cor- 
rections to total magnitudes as stated in Papovich et al. (2010). 
We matched the IRAC sources to the optical+IR data within 
1" radii. The final catalog covers a wavelength baseline of 
0.4-8 ^m. 

We use the photometric redshift catalog from Papovich 
et al. (2010) derived using EAZY vl .0 Brammer et al. (2008). 
In summary, we used the default galaxy spectral energy dis- 
tribution templates with the default A"-band prior based on the 
luminosity functions of galaxies in a semi-analytic simlua- 
tion. In addition to photometric redshifts, we measure the full 
photometric -redshift probability distribution function, P(z), 
normalized such that J P{z) dz = 1 when integrated over all 
redshifts. As discussed in Papovich et al. (2010), compared to 
existing spectroscopic redshifts, the best-fit photometric red- 
shifts have uncertainties of A([z sp - z p h]/[i +Z sp i) = 0.04 de- 
rived from the normalized median absolute deviation (Beers 
et al. 1990) in the range 1 < z sp < 2 (which includes the 
redshift-range of interest in this paper). 

We use stellar masses, star-formation rates (SFRs) and spe- 
cific SFRs (sSFR, the SFR per unit stellar mass), and rest- 
frame U - B and U — V colors as derived in Papovich et al. 
(2012). To summarize, we fit the 10-band galaxy photometry 
with model spectral energy distributions with Chabrier initial 
mass function and solar metallicity using stellar population 
synthesis models from Bruzual & Chariot (2003), allowing 
for a range of stellar population age, star-formation history 
and dust attenuation using the Calzetti et al. (2000) law. The 
rest-frame U-B and U — V colors and SFRs are derived us- 
ing the best-fit models, where we measure the SFR averaged 
over the prior 100 Myr. We found that the rest-frame colors 
for galaxies in our sample agree with those independently de- 
rived from the best-fitting templates for the photometric red- 



shifts using EAZY with a median rest-frame color differences 
of A(U-B) = 0.03 and A(U-V) = 0.04 mag with standard 
deviations a(A[U-B]) = 0.04 and a(A[U-V]) = 0.08 mag. 

We adopt for our galaxy sample the limiting magnitude 
from the Williams et al. (2009) catalog of K = 22.9 mag. 
This equates approximately to a 10cr detection limit. This 
magnitude limit for a galaxy at z = 1 .6 corresponds to a stel- 
lar mass of « 2 x 10 10 M , for a solar metallicity stellar 
population with a maximum mass-to-light ratio (formed in 
a burst at Zf = 5). Galaxies with younger stellar popula- 
tions will have lower mass-to-light ratios. Therefore, while 
our sample is "complete" for all galaxies with stellar masses 
above ps 2 x 10 10 M , our sample is sensitive to galaxies 
with younger stellar populations with stellar masses below 
this limit. Therefore, our stellar mass limit is "conservative" 
in the sense that we are sensitive to all galaxies above this 
mass limit. The choice for this limiting stellar mass is also 
motivated by our ability to recover accurate morphological 
information from the //S77CANDELS images with GALFIT. 
Galaxies at the stellar mass limit and with the fiducial stellar 
population above have magnitudes of m(F125W) = 24 mag, 
up to which we can recover robust quantitative morphological 
parameters with accuracies better than 40% (see below, § 3.1 
and the Appendix). We therefore adopt a limiting mass of 
2 x 10 10 Mq where we are complete for all galaxies, and for 
which our analysis is robust. 

We focus on the subset of the UDS with CANDELS imag- 
ing. In addition to the portion of the UDS covered by CAN- 
DELS, we find it is useful to include the full UDS dataset 
when making comparisons not involving size or morphology. 
This larger data set covers an area of 0.70 deg 2 (as discussed 
in Williams et al. (2009) and Papovich et al. (2012)) and is 
nearly 10 x larger in area than the CANDELS field. 

As in Papovich et al. (2012) we select galaxies around the 
redshift z = 1 .6 by defining the integrated photometric redshift 
probability distribution function, V z = J P(z) dz integrated of 
the redshift range z = z cen ± & with z cen = 1 .625 (Papovich 
et al. 2010, which corresponds to the mean redshift of the 
cluster in this field) and 5z = 0.05 x (1 +z C en), which is typi- 
cal of the photometric redshift 68% confidence limits at this 
redshift. Here we select all galaxies with V- = 0.4. As shown 
in Papovich et al. (2010), this V z limit is required to include 
blue, star-forming galaxies with spectroscopically confirmed 
redshifts at z ~ 1 .6, and it includes the red, quiescent galaxies 
at this redshift. Using a higher V z limit restricts the sample 
more to quiescent galaxies (as was used by Papovich et al. 
2012). Selecting galaxies this way has the advantage that it 
uses all information about the galaxies' photometric redshift 
probability distribution functions. In contrast, using samples 
of galaxies selected solely with a best-fit photometric redshift 
excludes information, and galaxies selected in some redshift 
interval will be very dependent on the errors in the photomet- 
ric redshifts and subject to biases. Regardless, in the case here 
we compare the properties of both cluster and field galaxy 
samples constructed from the same photometric -redshift se- 
lection and using the same CANDELS data set. Therefore, 
any systematics or biases resulting from the selection affect 
both samples equally. As a result, the relative comparison be- 
tween the galaxies in the cluster and field is robust. 

We select all galaxies with V- > 0.4. This selection pro- 
vides a sample of 433 galaxies with coverage in CANDELS 
imaging with mean (best-fit photometric) redshift z= 1.61 
with a standard deviation of 0. 12. 
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3. ANALYSIS 

3.1. Quantifying Galaxy Morphologies 

We used GALFIT (Peng et al. 2002) to fit parametric mod- 
els to the surface brightness profiles of our sample galaxies 
in the WFC3 F125W imaging. We chose the F125W band- 
pass for this analysis because it corresponds approximately 
to the fi-band in the rest-frame at z = 1.6. Many studies of 
the morphologies of distant galaxies (especially in clusters) 
are based in the rest-frame B-band (see, e.g., Blakeslee et al. 
2003; Homeier et al. 2005; Postman et al. 2005; Holden et al. 
2005; Blakeslee et al. 2006; Mei et al. 2009; Santos et al. 
2009) , and our choice of the F125W bandpass here allows for 
direct comparisons with these other studies. In addition, the 
CANDELS imaging allows for studies of galaxy morphology 
in the rest-frame B-band out to z ~ 3 using the F160W band- 
pass, which would permit the extension of our work to higher 
redshift. Therefore our choice here seems prudent. Regard- 
less, our tests (see below) show that none of our conclusions 
would be changed if we used the F160W bandpass instead. 

GALFIT models the surface brightness profile of a galaxy 
as exp(-B//? e ff) 1 /" (Sersic 1968), where « = 1 corresponds to 
a exponential (disk) profile and « = 4 corresponds to the de 
Vaucouleur's profile. GALFIT convolves these models with 
a user-defined PSF We used a PSF constructed using Tiny- 
Tim v7.2 (Krist 1995) and combined with the same dither po- 
sitions as the CANDELS data. For each galaxy we fit with 
GALFIT the effective semimajor axis (a e ffX background, axis 
ratio (b/a), total magnitude, and Sersic indices (n) as free pa- 
rameters. 

We compared our Sersic index measures using the F125W 
imaging for the galaxies in our z ~ 1.6 sample against in- 
dependent measurements derived using GALAPAGOS (Bar- 
den et al. 2012) with the CANDELS WFC3/F160W imag- 
ing as described by Bell et al. (2012) and van der Wei et al. 
(2012). We find a median offset in Sersic index of -0.03 dex 
in A(«) / n with a scatter of 0. 12 dex and we find a median off- 
set in effective radius of 0.01 in A(r e ff)/7 e ff with a scatter of 
0.06 dex. The uncertainties are nearly independent of galaxy 
magnitude, with only minor improvement for J < 23 mag 
compared to 23 < J < 24 mag. Because these differences are 
small, our conclusions would be unchanged if we used fits 
from the F160W imaging instead of the F125W imaging. 

We performed a series of simulations, inserting model 
galaxies into the HST data, and using GALFIT as described 
above to recover their parameters. Our analysis of these sim- 
ulations showed that the recovered effective radii are accu- 
rate to better than 40% and the Sersic indices to better than 
25% for galaxies with r e g = 2 kpc and n = 4 with magnitude 
m(F!25W) = 24 mag, near the stellar-mass limit of our sam- 
ple, (2 x 10 10 Mq, see § 2). Galaxies with parameters such 
as these represent the worst case scenario, and errors for mo st 
objects will be significantly smaller than this. The uncertain- 
ties are correlated between parameters, similar to the findings 
of Haussler et al. (2007). For more information, see the ap- 
pendix. 

3.2. Quantifying Galaxy Environment 

We tested different ways to characterize the environment 
of each galaxy in our z = 1.6 sample. We found that using 
samples based on the projected distance from the cluster cen- 
ter provided the most significant results. For this reason, we 
defined two samples. We selected a sample of "cluster" galax- 
ies with V z > 0.4 and projected distances ^ pro j < 1.5 Mpc 



of the cluster center (with the cluster center defined by Pa- 
povich et al. 2010). A subsample of "field" galaxies with 
^proj > 3 Mpc was selected for comparison. 

We also calculated environment measures based on the pro- 
jected distance to the Mh nearest neighbor, (e.g., Dressier 
1980). We computed the Mh nearest neighbor distance from 
each galaxy in our V z > 0.4, z ~ 1 .6 sample to (1) each other 
galaxy in the z ~ 1.6, V- > 0.4 sample, as well as to (2) 
each galaxy in a sample that have best-fit photometric red- 
shifts 1 .50 < z < 1 .75. Additionally, we tested the Mh nearest 
neighbor distance with the requirement that the stellar-mass 
ratio between each galaxy and its neighbors be greater than 
or equal to a certain value over a range from 1:1 to 1:10 (as 
proposed by Haas et al. 2011). In practice, we find that us- 
ing a nearest-neighbor distance estimator, e.g., the N = 2nd- 
nearest neighbor with D2 < 0.25 Mpc, identified the majority 
of galaxies with V z > 0.4 within R pro j < 1 Mpc of the clus- 
ter. However, all of the estimators selected only a small 
fraction (< 20%) of quiescent galaxies with V : > 0.4 within 
1 Mpc < R <1.5 Mpc from the cluster: the Mh nearest neigh- 
bor distances of these galaxies are simply too large to place 
them in high local number density regions based on this statis- 
tic. Haas et al. (2011) showed that including a requirement 
on the stellar-mass ratio between neighboring galaxies in the 
density estimator improves the strength of the correlation be- 
tween halo mass and density. However, our tests showed that 
this offers no improvement. Indeed, by requiring neighbor- 
ing galaxies to be at least as massive as the galaxy in ques- 
tion, nearly all of the most massive galaxies in the cluster are 
shifted to "low" density regions because there are no other 
galaxies of comparable mass in the sample (and therefore the 
distance to the Mh nearest neighbor with comparable mass 
becomes very large). 

It has been shown that Mh nearest neighbor measurements 
are good probes of internal cluster properties and galaxy- 
group-sized halos while fixed aperture methods better exam- 
ine the cluster as a whole (Muldrew et al. 2012). For this 
reason, Woo et al. (2012) conclude that using the relative 
projected distance of a satellite from its halo center is better 
suited to studying processes associated with clustering on the 
largest scales than Mh-nearest neighbor statistics. Moreover, 
we are motivated by the results from other studies that have 
found that the environment affects galaxies as they become 
satellites of larger dark-matter halos, whereas the evolution 
of central galaxies depends primarily on the halo mass (see 
§1). Therefore, here we adopt the projected distance from the 
cluster center as our measure of environment for our study. 
We define our samples as the "cluster" (with R pm j < 1 .5 Mpc) 
and "field" (with R pm j > 3 Mpc). In addition, we will divide 
the "cluster" sample spatially into two subsamples: the cluster 
center with R pm j < 1 .0 Mpc and the cluster outskirts defined 
as the annulus with 1.0 Mpc < R pm j < 1.5 Mpc. 

4. THE RELATION BETWEEN ENVIRONMENT, MORPHOLOGY, 
COLOR AND MASS AT Z = 1 .6 

Figure 1 shows the relation between the U - B rest-frame 
color and stellar mass for galaxies in our sample in the cluster 
(Bproj < 1.5 Mpc) and the field (R plo j > 3 Mpc), divided into 
subsamples of disk-dominated galaxies (Sersic index n < 2) 
and bulge-dominated galaxies (Sersic index n > 2). 

In all panels a strong U-B color-mass, and color-mass- 
morphology relation is apparent. Galaxies with high stel- 
lar mass have red colors (similar to, Williams et al. 2009; 
Quadri et al. 2012). The Sersic index is correlated with mass 
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Figure 1. The rest-frame U — B color - stellar mass relation for galaxies at z. = 1 -6. Each panel shows a subset of galaxies in the field and cluster with high and 
low Sersic index. The left column shows the relation for galaxies in the field, S plo j > 3 Mpc from the cluster center. The right column shows the relation for the 
cluster, with R„ m j < 1.5 Mpc. The upper row shows the relation for the subsample with low sersic indices (n < 2), and the lower row shows the subsamples with 
high Sersic indices (n > 2). Red-colored points denote objects that satisfy the color selection for quiescent galaxies. The thick blue lines denote the estimated 
stellar-mass completeness for stellar population with a maximal mass-to-light ratio, which formed in burst at z.f = 5 with solar metallicity and subsequent passive 
evolution. Our analysis is complete for all galaxies with stellar masses > 2 X 10 10 Mq, where our analysis shows our quantitative morphological parameters are 
robust. 

the star-forming "sequence" (e.g., Daddi et al. 2007; Noeske 
et al. 2007). Therefore, although strictly speaking the color- 
selection criteria in equation 1 separates galaxies exhibiting 
highly suppressed star-formation from those that do not. For 
brevity, we refer to these two samples as "quiescent" and 
"star-forming" galaxies respectively. 

For much of the analysis here, we use a one-sided 
Wilcoxon-Mann-Whitney rank-sum (WMW) test (Mann & 
Whitney 1947) to place a significance on the likelihood that 
the cluster and field galaxies have similar properties and are 
drawn from the same parent population. A summary of these 
likelihoods is given as probabilities (p-values) in Table 1 . 

4.1. Do the U — B Colors Depend on Environment? 

We do observe a slight increase in the U — B color of galax- 
ies associated with the cluster compared to the field. This is 
illustrated in Figure 2, which shows the cumulative distribu- 
tion of U - B color for both samples using the data for the 
CANDELS subset. The WMW test gives that the CANDELS 
galaxies in the cluster have higher U — B color, with a WMW 
significance p = 0.069. Using the larger full UDS sample im- 
proves this significance, with p = 0.033 (~ 2a). 

Most of the difference in the U—B colors stems from differ- 
ences in the colors of the star-forming galaxies. The WMW 
test for the U — B colors finds no difference in the distribu- 
tions for the quiescent galaxies (p > 0.3 for all samples). The 
star-forming galaxies in contrast show slight evidence that the 



and color. Exactly half (82/164) of bulge-dominated (n > 2) 
galaxies in our sample lie on the red sequence, and the ma- 
jority (238 /269 = 88.5%) of disk-dominated galaxies (n < 2) 
reside in the blue cloud. Disk-dominated galaxies populate 
the full range of U—B color. This is similar to the findings 
of previous studies at this redshift, for example, Wuyts et al. 
(201 1) and Bell et al. (2012), and references therein. 

Figure 1 also shows the color-mass-morphology distribu- 
tions for galaxies with colors consistent with older, passively 
evolving stellar populations. We selected these "quiescent" 
galaxies using a variant of the UVJ selection (Wuyts et al. 
2007; Williams et al. 2009) based on the observed z—J and 
/— [3.6] colors as defined by Papovich et al. (2012), 



(z-J)ab > 1.3 mag 
(7-[3.6]) AB <2.1mag 

(z-7) AB >0.5 + 0.55(y-[3.6])AB. 



(1) 



This color-selection classification is efficient in isolating 
galaxies with very different specific SFRs. In Papovich et al. 
(2012) we showed that using this color selection, the major- 
ity (~90) of quiescent galaxies selected have specific SFRs 
< 10~ 2 Gyr _1 . Therefore, while their SFR may not be zero, 
it is highly "suppressed" (e.g., Kriek et al. 2006). In con- 
trast, galaxies that do not satisfy our definition for "quies- 
cent" have specific SFRs ~ 1 Gyr" 1 , at least one order of 
magnitude higher, which are characteristic of galaxies on 
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U-B colors of the galaxies in the cluster are redder, although 
the significance is weak with p = 0.07 (about 1.5cr). There 
is little difference in this /9-value for the comparison between 
the smaller CANDELS samples and the larger UDS samples. 

If this difference in the U-B colors of the star-forming 
galaxies is real, then it would provide evidence that the star- 
forming galaxies associated with the cluster have more dust 
extinction, which reddens the colors. This could imply that 
the SFRs we derive from the SED fitting may be lower limits, 
missing highly extincted star-formation, which is indicated by 
the redder U-B colors. 

For reasons discussed below, we tested if galaxies at larger 
clustercentric distances, defined by the annulus 1 Mpc < 
^proj < 1 -5 Mpc, have different U-B colors compared to those 
galaxies within R pm j < 1 Mpc from the cluster. However, 
the WMW test gives no measureable difference with p-values 
> 0.11 for all subsamples, which is not significant. There- 
fore, the cluster galaxies within R pm j < 1 Mpc have similar 
U-B colors as those galaxies within 1 Mpc < R pm j < 1 .5 Mpc. 
We repeated this test for only the subsample of galaxies with 
M = (2-9) x 10 10 M Q , but find no substantial differences. 

4.2. Do the specific SFRs of Galaxies Depend on 
Environment? 

Considering the combined samples (quiescent and star- 
forming galaxies), the WMW does not rule out the hypothesis 
that the specific SFR distributions are drawn from the same 
parent populations (all have p > 0.22), and we see no dif- 
ference in the specific SFR distributions for more restrictive 
samples of quiescent or star-forming galaxies. We also see 
no significant difference using either the smaller CANDELS 
sample or larger UDS sample. There is no apparent difference 
in the specific SFR distributions between the cluster and the 
field, nor is there any significant difference in the distributions 
for galaxies in the cluster core {R woi < 1 Mpc) and those in the 
^proj = 1-1.5 Mpc annulus. 

On the surface this result appears to contrast some stud- 
ies that find an increase in the SFR (or SFR surface density) 
of galaxies in the cores of clusters at this redshift (e.g., Tran 
et al. 2010; Kocevski et al. 2011). For example, Tran et al. 
(2010) found an elevated IR luminosity projected surface den- 
sity associated with this cluster compared to the field. We at- 
tribute this apparent discrepency to the fact that in our study 
we have not included IR data, and therefore our SFRs (and sS- 
FRs) may be underestimated. Indeed, this is consistent with 
the star-forming galaxies in the cluster having redder U-B rest 
frame colors from §4.1. As discussed above, this may im- 
ply these galaxies have high extinction, which would imply 
higher extinction-corrected SFRs that we are missing in our 
measurements. 

4.3. Do the Stellar Masses Depend on Environment? 

In general, there is no substantial difference in the distribu- 
tion of stellar masses when all galaxies (quiescent and star- 
forming) are considered in the field and the cluster. For the 
more restrictive sample of quiescent galaxies only, those as- 
sociated with the cluster show a possible increase in stellar 
mass relative to the field. The WMW test gives a signifi- 
cance p = 0.065 for the CANDELS samples. However, this in- 
creases to p = 0.078 when all galaxies in the larger UDS sam- 
ples are considered, implying a weaker significance. There- 
fore, there is slight evidence (1 .5a significance) that quiescent 
galaxies associated with the cluster have larger stellar masses 
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Figure 2. The cumulative distribution of rest-frame U—B colors for all 
galaxies in the field (fi pl0 j > 3 Mpc, solid blue line) and in the cluster 
(^pioj < 1-5 Mpc, dashed red line). Each curve is normalized to the total 
number of galaxies in that sample (these curves include all galaxies in the 
UDS, not only those in the CANDELS field). The left panel shows the 
distribution for all quiescent galaxies with stellar masses > 2 X 10 10 Mq. 
The circle points on each curve show the median value of the distributions. 
There is weak evidence that the galaxies associated with the cluster have 
a median U — B color that is A(t/ - B) ~ 0. 1 mag redder than those in the 
field. The WMW gives p = 0.033 for the total UDS samples, but declines 
to p = 0.069 using only the CANDELS samples. The right panel shows a 
subset with moderate stellar mass 2 X 10 10 -9 X 10 10 Mq. The WMW test 
finds a reduced significance, p > 0.100 (see Table 1). 

compared to those in the field. The difference appears driven 
by the higher stellar masses of the quiescent galaxies in the 
cluster compared to the field, which is evident in Figure 1 . 

The stellar masses of galaxies in the 1 Mpc < R pm j < 1.5 
Mpc annulus show evidence of being different from those 
with R pm j < 1 Mpc. Considering all galaxies (star-forming 
and quiescent), the galaxies within R pm j < 1 Mpc have stellar 
masses approximately a factor of 1 .9 higher than those within 
1 Mpc < R pmi < 1.5 Mpc (with WMW p = 0.0053). This sig- 
nificance is reduced (p = 0.078) when considering only qui- 
escent galaxies. Interestingly, there is no difference between 
the galaxies within 1 Mpc < R pro j < 1.5 Mpc and those in the 
field, R pm j > 3 Mpc, (p > 0.2), suggesting galaxies in this an- 
nulus share a mass distribution more characteristic of the field 
at this redshift than the highest-density region associated with 
the cluster. This is possible evidence that these galaxies were 
recently members of the field and are in the process of being 
accreted into the cluster environment. 

4.4. Do the Sersic Indices of Galaxies 
Depend on Environment? 

We now compare the difference in Sersic indices between 
the cluster and field populations in more detail. This is il- 
lustrated in Figures 3 and 4. Figure 3 shows the astrometric 
positions of quiescent and star-forming galaxies at z = 1 .6 in 
the CANDELS field. The symbols correspond to the Sersic 
index of the quiescent galaxies. Figure 4 shows the Sersic in- 
dices of galaxies in the CANDELS field as a function of their 
clustercentric distance. These figures show that there appears 
to be an excess of galaxies with n < 2.5 in the annulus 1 Mpc 
< ^proj <1.5 Mpc from the cluster relative to either the sam- 
ple with R pm j < 1 Mpc or the field sample (R pm j >3 Mpc), 
and that this effect is driven by intermediate mass galaxies. It 
is also apparent that quiescent galaxies with n < 1 are only 
found in the 1 Mpc < R pro j <1,5 Mpc annulus. As discussed 
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Table 1 

Summary of Wilcoxon-Mann- Whitney U probabilities 
comparing the properties of z= 1.6 cluster and field galaxies in the CANDELS UDS field 



Sample 


^cluster A/field 


p(U-B) 


p(logM„) 


p(log sSFR) 


pfeff) 


p(n) 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


All Galaxies 


logM*/M© > 10.30 












UDS 


56 2061 


0.0331 


0.125 


0.296 








37 110 


0.0692 


177 




422 


n ns on 


All Galaxies 


10.30 < log M*/M© 


< 10.95 










UDS 


43 1701 


0.100 


0.308 


0.412 






CANDELS 


27 94 


0.229 


0.390 


0.263 


0.493 


0.091 


Quiescent Galaxies 


logM»/M Q > 10.30 












UDS 


35 1083 


0.458 


0.0781 


0.216 






CANDELS 


21 51 


0.397 


0.0646 


0.202 


0.0571 


0.0234 


Quiescent Galaxies 


10.30 < logM«/M < 10.95 










UDS 


24 870 


0.291 


0.365 


0.292 






CANDELS 


13 41 


0.292 


0.310 


0.420 


0.107 


0.0164 


star-forming galaxies 


logM»/M Q > 10.30 












UDS 


21 978 


0.0699 


0.466 


0.409 






CANDELS 


16 59 


0.0737 


0.339 


0.193 


0.280 


0.238 


star-forming galaxies 


10.30 < logM./Mg 


< 10.95 










UDS 


19 831 


0.0586 


0.353 


0.444 






CANDELS 


14 53 


0.137 


0.254 


0.297 


0.181 


0.457 



Note. — ( 1 ) Description of samples used, (2) number of cluster galaxies in WMW test, (3) number of field galaxies used in WMW 
test. Other columns are the WMW significance of the difference in the cluster and field samples between the distributions of (4) stellar- 
mass, (5) specific SFR, (6) effective radii, (7) Sersic indices. Low probabilities (p-value) indicate more significant differences in the 
distributions. 



in § 3.1, to test if this is a result of our choice of F125W band- 
pass or analysis, we checked our Sersic indices against those 
computed independently from the F160W data from Bell et al. 
(2012), and we found no change. 

4.4. 1 . Comparison Between the Cluster and the Field 

Table 1 shows the WMW p-values comparing the proper- 
ties of galaxies in the cluster and field environments. Based 
on our analysis there is no evidence that star-forming galaxies 
in the field and cluster have differences in their Sersic indices. 

We do observe differences in the distributions of Sersic in- 
dices for quiescent galaxies. There is an envelope in the Ser- 
sic index and stellar mass distribution for quiescent galaxies: 
more massive quiescent galaxies generally have higher Sersic 
indices, and only a few of the most massive (> 9 x 10 10 M@) 
galaxies have n < 2 (Figure 1). Most, 81/108 = 75%, of the 
quiescent galaxies in both the field and cluster have high Ser- 
sic indices, n > 2, and are bulge-dominated (consistent with 
Bell et al. 2012; Wuyts et al. 2011; Papovich et al. 2012). 
However, more notably, the fraction of galaxies with low Ser- 
sic index « < 2 that are quiescent is significantly higher in 
the cluster (6/16=38%) compared to the field (8/43 = 15% ), 
for galaxies above our mass limit. This is significant at about 
2.5tr from binomial statistics, although it is based on small 
numbers of galaxies. As illustrated in Figures 3 and 4 and as 
we discuss below, most of this difference comes from galax- 
ies in the annulus 1 Mpc < R WO] < 1 .5 Mpc from the cluster, 
and suggests some mechanism suppresses star-formation in 
galaxies with low Sersic index at large clustercentric radii. 

The WMW probabilities in Table 1 show the evidence for 
differences in the Sersic indices of galaxies associated with 
the cluster and the field. Most of this difference is isolated 
to the quiescent galaxies: star-forming galaxies in both the 
field and cluster have median Sersic indices n me( j ~ 1 .4, with 
no significant differences, p > 0.2 . Quiescent galaxies as- 



sociated with the cluster have lower median Sersic indices, 
«med = 2.6 versus n me d = 3.3 in the field, and the WMW test 
shows this is significant with p = 0.023 (approx. 2a) consid- 
ering the full sample with M > 2 x 10 10 M & . Figure 5 shows 
the cumulative distribution of the Sersic indices of quiescent 
galaxies in the field and the cluster. 

The right panel of Figure 3 shows the cumulative distribu- 
tions for the more restrictive subsamples of quiescent galaxies 
with moderate stellar mass, 2 x 10 10 -9 x 10 10 M Q . Here the 
median Sersic index is « me d = 2.4 for the cluster galaxies, and 
"med = 2.8 for the field galaxies, and the significance increases 
to p = 0.016. 

This result has a significance of « 2er, and is based on rela- 
tively small samples: 21 quiescent galaxies with stellar mass 
> 2 x 10 10 M Q which lie within R pm j < 1 .5 Mpc from the clus- 
ter as covered by CANDELS (table 1). To test further the sig- 
nificance of this result, we used a bootstrap Monte Carlo test 
to estimate a likelihood that we obtained the observed differ- 
ence in Sersic indices by chance. We used as a parent sample 
the combined samples of all quiescent galaxies in the cluster 
and field. We then constructed 10,000 pairs of cluster-sized 
and field-sized samples of A^(cluster) and A^(field) objects ran- 
domly selecting from the parent sample (with replacement), 
where A^(cluster) and A^(field) are equal to the size of the orig- 
inal cluster and field samples, respectively. We then recom- 
puted the WMW p value between each pair of random cluster 
and field samples. In this test, we obtained a p value as small 
(or smaller) than the one we observe in only 0.6% of the ran- 
dom samples when drawing from the full sample of galax- 
ies with M > 2 x 10 10 M Q (this increases slightly to 0.9% of 
cases for M = 2 x 10 10 -9 x 10 10 M© moderate-mass sam- 
ples). Therefore, the likelihood that we would have obtained 
our result by chance if there were no actual difference be- 
tween the cluster and field objects is smaller than 1% (equiv- 
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Figure 3. The spatial distribution of galaxies at z = 1.6 in the CANDELS UDS field. Symbols denote quiescent (large symbols) and star-forming galaxies (small 
boxes), as labeled. For the quiescent galaxies, the symbol type corresponds to Sersic index, as indicated in plot legend. The concentric dashed circles show 
projected distances of R pro j =1, 1.5, and 3 Mpc from the astrometric center of the cluster using the coordinates of Papovich et al. (2010). The quiescent galaxies 
in the annulus 1 Mpc < S plo j < 1.5 Mpc have lower Sersic indices and larger effective sizes compared to other quiescent galaxies. 



alent to the probability of > 2.6<7 significance for a Gaussian 
distribution). 

4.4.2. Comparison Between the Cluster Core and the 
1 Mpc < R pm j < 1.5 Mpc Annulus 

As shown in Figure 5, the quiescent galaxies in the 1 
Mpc < R pm j < 1 .5 annulus have lower median Sersic indices 

n rned 5MpC = 2.1, compared to those within 1 Mpc of the center 
of the cluster, which have median «,^g^ Ipc = 3.0. Indeed, four 
out of the six moderate mass quiescent galaxies in the annu- 
lus have Sersic indices n < 2.5, including three galaxies with 
n < 1 . This is evident in Figures 3 and 4, which show that 
these low-Sersic-index quiescent galaxies are found preferen- 
tially in at 1-1.5 Mpc from the cluster. This difference is even 
more apparent considering the most restrictive moderate-mass 
subsample with 2 x 10 10 -9 x 10 10 M Q . For this sample, the 
median Sersic indices are n= 1.0 and n = 2.4 for the annulus 
and cluster core respectively. 

The significance of this result is somewhat limited by the 
small number of objects in the samples: there are only 95 
quiescent galaxies in our cluster and field samples combined. 
To measure the significance of this result, we considered two 
Monte Carlo bootstrap simulations to determine how often we 
would obtain this result by chance. 

We first calculate how often we would obtain the results by 
chance if the quiescent galaxies in the 1-1.5 Mpc annulus are 



randomly drawn from the parent field sample. For this test we 
constructed two random samples, one having N(l-1.5Mpc) 
members (the number of quiescent galaxies in the 1-1.5 Mpc 
annulus) and the other having N(field) members. Both sam- 
ples are drawn from the quiescent field population with re- 
placement. From this bootstrap test, we found that the WMW 
p value was smaller than the measured value in only 0.44% of 
the cases, which corresponds to a significance of « 2.6er (as- 
suming a Gaussian distribution). Therefore, it seems that even 
with the small sample sizes here, the lower Sersic indices of 
the quiescent galaxies in the annulus 1 Mpc < R WOj < 1 .5 Mpc 
are statistically significant compared to quiescent galaxies in 
the field. 

We next repeated this bootstrap test to estimate how often 
we would obtain the results by chance if the quiescent galax- 
ies in the annulus are randomly drawn from the parent cluster 
quiescent sample. Again, we constructed random samples, 
one having N(l-1.5 Mpc) members, and one having N(<1 
Mpc) members. Here, we find that the likelihood that the qui- 
escent galaxies in this annulus are drawn from the same pop- 
ulation as those in the rest of the cluster is p\, { 



iootstrap 



= 0.039 

(« 1.8er) primarily because of the smaller number of quies- 
cent galaxies in the cluster sample. Larger samples of quies- 
cent galaxies (especially those in clusters) at this redshift are 
required to improve the statistical significance of this obser- 
vation. 
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To summarize this section, there is evidence that the qui- 
escent galaxies in the cluster have lower average Sersic in- 
dices compared to quiescent galaxies in the field. Further- 
more, there is evidence that this result is driven by a popula- 
tion of quiescent galaxies in an annulus 1-1.5 Mpc from the 
cluster. Members of this population are found to have Sersic 
indices inconsistent with random samples of quiescent galax- 
ies taken from the cluster core (R pm i < 1 Mpc) or from the 
field. Instead, they are morphologically more similar to star- 
forming galaxies found in all environments. 

4.5. Do the Effective Radii of Galaxies 
Depend on Environment? 

Because the quiescent galaxies in the cluster show evi- 
dence for smaller Sersic indices, we also compared the ef- 
fective radii of these galaxies to those in the field. The me- 
dian effective radius of quiescent galaxies with stellar mass 
> 2 x 10 10 M associated with the cluster is rf^ s ^ A = 2.0 kpc 
compared to those in the field which have f^ff^ed = 1-4 kpc. 
(We remind the reader that these radii are the circularized ef- 
fective radius). The WMW test gives a probability, p = 0.057. 
This is nearly identical to the findings of Papovich et al. 



(2012), who used the same data, but with slightly different 
selection criteria and mass limits to focus on the quiescent 
galaxies (see also § 2). Figure 6 shows the cumulative distri- 
bution of the effective radii of quiescent galaxies in the field 
and cluster, which illustrates that the cluster galaxy sizes are 
larger on average. 

We also considered the effective sizes of the quiescent 
galaxies within R pm j < 1 Mpc of the cluster center compared 
to those in the annulus 1 Mpc < R pm j < 1.5 Mpc because 
these show evidence for having different Sersic indices. Fig- 
ure 6 shows the cumulative distributions of these subsamples. 
The quiescent galaxies in the 1 Mpc < R pro j < 1.5 Mpc annu- 



lus have median effective radii, r, 



1-1.5 Mpc 

eff.med 



2.7 kpc, about 1.3 



times larger than other quiescent galaxies in the cluster or the 
field. 

Because of the relatively small size of our galaxy samples, 
we again used a Monte Carlo bootstrap simulation (similar to 
the previously described simulations) to test the significance 
of these results. Based on this test, the likelihood that the 
quiescent galaxies in this annulus are drawn from the cluster 
sample is 9.6%. Similarly, the likelihood that the quiescent 
galaxies in this annulus are drawn from the same population 
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The results of this section can be summarized as fol- 
lows: we find that quiescent galaxies associated with the 
cluster have low Sersic indices and possibly larger effective 
radii (when compared with quiescent galaxies in the field) 
consistent with disk-like morphologies (« ~ 1). This effect 
is driven by galaxies located in an annulus 1-1.5 Mpc from 
the center of the cluster, whereas quiescent galaxies in the 
cluster core and field exhibit Sersic indices more typical 
of a bulge dominated morphology. Furthermore, the most 
massive galaxies are located in the cluster core while stellar 
masses of galaxies in the annulus are comparable to those 
found in the field. This suggests that the sample of quiescent 
galaxies found 1-1.5 Mpc from the center of the cluster is 
"contaminated" with recently quenched star-forming galaxies 
accreted from the field. 



Figure 5. The cumulative Sersic-index distribution for quiescent galaxies in 
the field (solid blue line) and in the cluster (dashed red line). The curves 
for the total cluster and total field samples are normalized to the number 
of galaxies in each sample. The curves for the cluster subsamples ("field" 
and "annulus") are normalized by the number of galaxies in the total cluster 
sample. The left panel shows the distribution for all quiescent galaxies with 
stellar masses > 2 X 10 10 Mq. The left panel shows a subset with moderate 
stellar mass 2 X 10 10 — 9 X 10 10 Mq. The circular points on each curve show 
the median value of the distributions. The quiescent galaxies associated with 
the cluster have smaller median Sersic indices, and this is most pronounced 
in the moderate mass subsample. The green (purple) dashed curve shows 
the contributions to the cluster cumulative distribution from quiescent clus- 
ter galaxies in the cluster "center" with R < 1 Mpc ("annulus" with 1 Mpc 
< R <1.5 Mpc). In both panels, the lower Sersic indices of quiescent cluster 
galaxies is driven by the quiescent galaxies associated with the cluster in an 
annulus 1 Mpc < R <1.5 Mpc, which have the lower median values. 
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Figure 6. The cumulative distribution of effective radii for quiescent galaxies 
in the field (i? pro j > 3 Mpc, solid blue line) and in the cluster (ifprqj < 1.5 Mpc, 
dashed red line). Each curve is normalized to the total number of galaxies in 
each sample. The left panel shows the distribution for all quiescent galaxies 
with stellar masses > 2 X 10 10 Mq. The right panel shows a subset with 
moderate stellar mass 2 X 10 10 -9 X 10 10 Mq. The quiescent galaxies asso- 
ciated with the cluster generally have larger median effective radii. In each 
panel the green (purple) dashed curve shows the contribution to the cluster 
cumulative distribution from quiescent galaxies in the cluster "center" with 
R < 1 Mpc ("annulus" with 1 Mpc < R <1.5 Mpc). The circle points show 
the median value for each distribution. Quiescent galaxies associated with 
the cluster in an annulus 1 Mpc < R <1.5 Mpc have larger median effective 
radii by factors of 2 compared to the field. 

as those in the field is 8.3%. Given the relatively small 
sample size, it seems that larger samples would be able to 
further test these trends. 



5. DISCUSSION 

5.1. The Relation Between Color, Stellar Mass, and 
Morphology in the Cluster and Field at z= 1.6 

We study how the relationship between galaxy color, stellar 
mass, and morphology depend on galaxy environment by ex- 
amining these relations in the field and in a forming cluster at 
z ~ 1.6. We find evidence of a possible correlation between 
galaxy density with color and mass for galaxies with stellar 
masses > 2 x 10 10 M Q : compared to galaxies in the field, 
the cluster galaxies (with R pm j < 1 .5 Mpc) have redder rest 
frame U-B colors and higher stellar masses. The color differ- 
ence appears driven by star-forming galaxies (we find no dif- 
ference in the rest-frame U-B colors of quiescent galaxies). 
We interpret this as evidence for increased dust-obscured star- 
formation in galaxies closer to the cluster, supported by the 
higher IR-luminosity density in these galaxies found by Tran 
et al. (2010). The stellar mass difference appears driven by 
the more massive quiescent galaxies with > 9 x 10 10 M Q (see 
§ 4.2 and Figure 4), and we interpret this as evidence that the 
most massive galaxies at this redshift in this field are associ- 
ated with the cluster. 

There are morphological differences in the quiescent galaxy 
populations associated with the cluster and the field. Combin- 
ing our observations with the results from our previous study 
(Papovich et al. 2012), we find two main differences between 
the morphologies of quiescent galaxies in the field and that of 
the cluster. First, Papovich et al. (2012) found a lack of "com- 
pact" quiescent galaxies with r e s ~ 1 kpc in the cluster while 
such objects are present in the field. Similarly, many of these 
quiescent galaxies in the cluster show evidence for extended 
disks. Our analysis here reproduces the findings of Papovich 
et al. (2012), and here we report evidence that the Sersic in- 
dices of the quiescent galaxies in the cluster are lower than for 
such galaxies in the field. 

5.2. The Nature of Quiescent Galaxies in the Annulus 
flproj = 1-1.5 Mpc 

The population of quiescent galaxies in the annulus defined 
by projected distances 1 Mpc < R pm j < 1.5 Mpc from the 
cluster show the most differences in morphology compared 
to the field. As discussed in § 4.4, galaxies in this annulus 
have higher r e ff and lower Sersic indices. Indeed, consider- 
ing the full subsample with M > 2 x 10 9 M Q , their median 
values are ri_i.5M pc = 2.7 kpc and «i-i.5Mpc = 2.1 (compared 
to medians r<i Mpc = 1-5 kpc and n<i Mpc = 3.1 in the cluster 
center). The quiescent galaxies in the /? prc .j= 1-1.5 Mpc an- 
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Figure 7. Examples of quiescent galaxies in the annulus 1 Mpc < R plo j < 
1.5 Mpc and low Sersic index (n < 2) and other galaxy subsamples. Each 
panel spans 12" on a side (about 100 kpc at z=1.6), and the colors show the 
HST F125W, (F125W+F160W)/2, and F160W as blue, green, and red, re- 
spectively. These images illustrate that these quiescent galaxies have similar 
morphologies as star-forming galaxies in both the cluster and field selected 
with low Sersic indices and over a similar range of stellar mass. The top 
row shows quiescent galaxies in the core with similar stellar masses but more 
typical, higher Sersic indices. 



nulus are more typical of star-forming galaxies in the field and 
cluster. Figure 7 shows examples of these quiescent galax- 
ies in the cluster annulus with low Sersic index (n < 2), and 
compares these with examples of star-forming galaxies with 
similar masses and Sersic indices (n < 2). These examples 
illustrate that the quiescent galaxies have morphologies simi- 
lar to some star-forming galaxies. In contrast, the figure also 
shows quiescent galaxies in the cluster core with more typi- 
cal, higher Sersic indices (n > 2). These have very different 
morphologies when compared to annulus quiescent galaxies 
with low Sersic indices. The star-forming galaxies have an 
interquartile range of 2.2 - 4.7 kpc in effective radius and 0.8 
- 3.3 in Sersic index, consistent with morphologies of these 
quiescent galaxies in the annulus. Because these galaxies ap- 
pear morphologically distinct, we speculate about their origin 
and fate. 

Bell et al. (2012) have argued that the Sersic index corre- 
lates the strongest with increasing rest-frame U — V color. We 
reproduce this relation in Figure 8, which shows the Sersic 
index versus U — V color for all z ~ 1 .6 galaxies in our CAN- 
DELS data with stellar mass > 2 x 10 10 M Q . 

We find a similar result to that of Bell et al., as most quies- 
cent galaxies lie above n > 2. However, many of the quiescent 
galaxies associated with the 1 Mpc < R pm j < 1.5 Mpc annulus 
appear as outliers. A high fraction («50%) of the quiescent 
galaxies in the annulus 1 Mpc < R pro j < 1 .5 Mpc have n < 2, 
compared to <15% of quiescent galaxies overall. This leads 
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Figure 8. The Sersic indices, n, versus the rest-frame U—V color for galaxies 
with stellar mass > 2 X 10 10 Mq. Bell et al. (2012) have argued that Sersic 
index best correlates with increasing U — V color. We reproduce that result 
here, as most quiescent galaxies lie above n > 2, indicated by the dashed line. 
Quiescent and star-forming galaxies in the subsample with 1 Mpc < R <1.5 
Mpc of the cluster center are indicated with large circles. A high fraction 
(Ri50%) of the quiescent galaxies in this annulus have n < 2, compared to 
<15% of quiescent galaxies overall. 

to the differences in the Sersic index distributions we see for 
quiescent galaxies in this annulus compared to the other sam- 
ples (see § 4.4). 

5.3. Environmental Effects on Galaxy Evolution at z = 1 .6 

At z = 1.6 there appear to be two evolutionary channels 
for quiescent galaxies, one that depends on the mass of the 
galaxy, and another that depends on environment. Mass is 
likely related to one channel of quiescence because both the 
field and cluster samples exhibit passive, bulge dominated 
galaxies (« > 2) which show no differences in their specific 
SFRs, or rest-frame U—B colors (see Table 1). This mass re- 
lated channel therefore has only a weak dependence on envi- 
ronment (cluster versus field) and is likely driven by processes 
related to the galaxy halo mass. This is similar to the findings 
of Peng etal. (2010). 

The environment also has some role in the morphological 
evolution of galaxies. One difference between the cluster 
and field quiescent galaxies is in their sizes. Cluster galax- 
ies have slightly larger effective radii than quiescent galax- 
ies in the field (see also Papovich et al. 2012) significant at 
~ 1.5a. There are also indications that the cluster galaxies 
have a higher "dry" merger rate (Lotz et al. 2013), which has 
been interpreted as evidence that the cluster galaxies have an 
increased or accelerated evolution. This view is consistent 
with the evolution of sizes of cluster galaxies, and the clus- 
ter luminosity function (Papovich et al. 2012; Rudnick et al. 
2012). The higher density region associated with the cluster 
appears to increase the merger rate of quiescent galaxies over 
that expected solely in the field. In the case of minor mergers, 
this has the effect of increasing their radii with only modestly 
increasing their mass. This process appears to work within 
r? pro j < 1 — 1.5 Mpc from the cluster center. 

The environment also seems to affect quiescent galaxies at 
higher projected distances from the cluster (1 Mpc < R plo j < 
1.5 Mpc). It is in this annulus where we see indications that 
the quiescent galaxies have larger effective radii and lower 
Sersic indices (see § 4.4 and Figure 4). The mass distribution 
of the quiescent galaxies in this annulus is also different from 
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quiescent galaxies in the cluster center, which suggests some 
process affecting galaxies accreted by the cluster is linked to 
the formation of these galaxies. 

As discussed above, the morphologies of the quiescent 
galaxies in this annulus are more typical of the morpholo- 
gies of star-forming, disk-dominated galaxies. Interestingly, 
the projected distances of the galaxies in this annulus lie out- 
side the upper limits of the virial radius estimated from galaxy 
velocity dispersions (~ 0.9 Mpc, under the assumption the 
galaxies are fully viralized, Papovich et al. 2010) and analysis 
of the X-ray data (Pierre et al. 201 1, ~ 0.5 Mpc). If our find- 
ings are confirmed by larger samples to be a general property 
of cluster galaxies at this epoch, it is implied that there is a 
mechanism which induces quiescence at large distances from 
the cluster center but leaves galaxy morphologies unchanged. 

One explanation for these observations is that these disk- 
dominated quiescent galaxies at projected distances 1 Mpc 
< ^pioj < 1-5 Mpc have been preprocessed in smaller group- 
sized halos, which are now merging with the cluster (e.g., 
McGee et al. 2009). Figure 3 shows that the late-type qui- 
escent galaxies may be located in ~2 groupings. Moreover, 
qualitatively many of the 1 < n < 2.5 quiescent galaxies in 
the field are not randomly distributed, but rather exist in small 
groups which are identified using Nth nearest neighbor mea- 
sures. 1 

However the idea of preprocessing does not explain why 
the morphological properties of the quiescent galaxies at 1 
Mpc < Rproj < 1 .5 Mpc would be significantly different from 
quiescent field galaxies. We favor an explanation where qui- 
escent galaxies at projected distances of 1 Mpc < R pm j < 1.5 
Mpc began as star-forming, disk-dominated galaxies before 
they entered the proximity of the forming galaxy cluster. At 
these distances they undergo truncation of their star forma- 
tion without disrupting their morphological profiles or sizes. 
At later times they will enter the cluster "halo", and, because 
they are quiescent, likely will undergo dissipationless merg- 
ers with other quiescent galaxies. This will then reconfigure 
their morphologies, consistent with the enhanced merger rate 
of Lotz et al. (2013). The mechanism for this transition from 
star-forming disk galaxy to quiescent galaxy appears to be- 
gin outside the cluster virial radius. If this occurs in galaxy 
groups, then it affects the galaxies only as they approach the 
region of the cluster. 

This interpretation may be explained by the expectation 
from theory that as gas falls into the region of large halos it 
is shocked to the virial temperature of the halo (e.g., White & 
Rees 1978; Birnboim & Dekel 2003). Recent cosmological 

1 Characterizing overdensity in our sample using the Mh nearest-neighbor 
method, we measured no appreciable differences in galaxy properties, using 
either the z. = 1 -6 sample defined by V- > 0.4 or the sample selected using the 
best-fit photometric redshifts in the range 1.5 < z < 1.75, This may be due 
to the fact that our samples rely largely on photometric redshifts, which is 
necessary to be complete for red galaxies at z ~ 1.6. Samples based on pho- 
tometric redshifts will suffer from projection effects along the line of sight 
(even the "good" photometric redshifts of our data, Az/(1 +z) = 0.04, yield 
line-of-sight proper distance uncertainties of 27 Mpc). This will smear any 
real correlations between galaxy density (environment) and galaxy properties. 
This is mitigated in lower redshift samples based on spectroscopic redshifts 
(e.g., Kauffmann et al. 2004; van der Wei et al. 2008; Peng et al. 2010; Cooper 
et al. 201 1). We speculate that the reason we are able to detect differences in 
the galaxy properties comparing the cluster to the field is that the overdensity 
of the cluster overcomes the smearing along the line-of-sight from photo- 
metric redshifts. Nevertheless, this smearing will affect our results, reducing 
any intrinsic trends between galaxy properties and environment. Therefore, 
it also follows that the results here show the minimum effect at z = 1.6 of en- 
vironment on galaxy properties such as color, masses, SFRs, specific SFRs, 
sizes, and Sersic index. 



Af-body and hydrodynamical simulations show that the veloc- 
ity shock from the virialized region extends beyond the virial 
radius (Cuesta et al. 2008; Birnboim et al. 2007; Dekel et al. 
2009), and this has effects on both the baryonic gas fraction, 
and the velocities and temperatures of the "shocked" infalling 
gas out to 1.5- 2R v j r as well as the stripping of subhalos at 
3-4i? vir . (Kravtsov et al. 2005; Bahe et al. 2012, E. Zinger, 
A. Dekel, et al., in preparation; P. Behroozi, R. Wechsler, et 
al. in preparation). 

The quiescent disk-dominated galaxies we observe at 1 - 
1.5 Mpc from the cluster may be inside this hot (shocked) 
medium, where gas accretion to them is then shut off (Dekel 
& Birnboim 2006; Croton et al. 2006). This process is of- 
ten referred to as gas "strangulation" (Balogh et al. 2000). 
It is possible these galaxies are on their first infall into the 
cluster: the cluster-crossing time is 5 Gyr (assuming the 1 
Mpc radius and typical peculiar velocity of «400 km s -1 , 
Pierre et al. 2012), and this exceeds the age of the Universe at 
z = 1.6. Under the assumption that the galaxies are on their 
first pass through the cluster (and have not already passed 
through), effects from a virial-shock-like process seem fa- 
vored because the galaxies lie outside the central cluster re- 
gion. Other processes that could affect these galaxies such 
as ram-pressure stripping (Gunn & Gott 1972; Abadi et al. 
1999), tidal stripping (Read et al. 2006), and galaxy harass- 
ment (Moore et al. 1996) are expected to be more significant 
nearer the central cluster regions where the galaxy and gas 
densities are higher. Because it is less likely that the quies- 
cent disk-dominated galaxies at 1-1.5 Mpc have had suffi- 
cient time to pass through the cluster, it seems more likely 
that they have consumed their gas supply and gas accretion 
has been "strangulated" by the hot medium. 

We also cannot rule out the possibility, however, that these 
galaxies have already passed through the outskirts of the 
larger dark matter potential and are now located at the turn 
around radius of the cluster. By passing through the cen- 
tral halo at large radii, the baryonic content of this popula- 
tion could have been altered without significantly restructur- 
ing their morphologies. Although we estimate that the cross- 
ing times for this cluster may be longer than the age of the 
universe at this redshift, this structure has grown since earlier 
epochs, and as it does not yet appear virialized the interpreta- 
tion of the kinematics is unclear. 

Another mechanism that may explain this processing of 
galaxies at large clustercentric distances is the assembly bias 
of dark matter halos (Gao et al. 2005; Wechsler et al. 2006; 
Croton et al. 2007). Simulations show that halos with early 
formation redshifts are more strongly clustered than halos of 
comparable mass, which form later. Assembly bias in the 
nearby universe (0.01 < z < 0.3) has been studied in the SDSS 
by Cooper et al. (2009) who find that galaxies with older stel- 
lar populations and higher stellar metalicities are preferen- 
tially found in higher density environments. This suggests 
that, at a given mass, early type galaxies in dense environ- 
ments were formed earlier, and thus ceased forming stars at 
an earlier epoch. This idea is supported by the comparison 
of ages of low redshift early-type galaxies between low and 
high local number density regions. Massive galaxies in high 
local number density regions appear to form their stellar pop- 
ulations in the range between z ~ 2 and z ~ 5, while their 
low local number density counterparts are 1-2 Gyr younger 
(Thomas et al. 2005). 

Simulations show that a cause for this assembly bias is that 
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Figure 9. The fraction of quiescent galaxies, /g, as a function of pro- 
jected distance from the center of the cluster. Galaxies have stellar mass 
> 2 X 10 10 Mq. Each datum shows the ratio of the number of quiescent 
galaxies to the total in each bin of projected distance. This fraction is printed 
above each datum. The horizontal dashed line shows the quiescent fraction of 
all galaxies in the same (including both cluster and field samples), /g = 0.52. 
The long-dashed curve shows the probability of obtaining the measured frac- 
tion or greater assuming a binomial distribution assuming an intrinsic fraction 
of /g. The dotted lines show significance levels of 1 and 2a. The binomial 
probability is very dependent on the intrinsic fraction of /g, and the value we 
assume is a conservative upper limit for reasons in the text. At projected dis- 
tances R pro j < 1 Mpc, we observe a significant enhancement of the quiescent 
fraction, but this does not extend beyond 1 Mpc. 

group-sized halos must compete for mass (both baryonic and 
dark matter) when in close proximity to much larger dark mat- 
ter halos. In effect, the group halos will be starved of dark 
matter and the galaxies they contain will be similarly starved 
of baryonic gas. These galaxies will exhaust their supply of 
star-forming gas more quickly than galaxies located in simi- 
lar halos in isolation. As this starvation is not directly asso- 
ciated with mergers or other more direct galaxy interactions, 
it is likely that the cessation of star formation will not dis- 
rupt a galaxy's late-type morphology. Our observations at 
z = 1 .6 may be due to this same affect: differences we see 
between field galaxies and those in the outskirts of the clus- 
ter may be due to an earlier formation of dark matter halos 
in regions of high density. Bell et al. (2012) studied galaxies 
over a large range of redshift in the CANDELS UDS field, 
finding the z ~ 1 .6 redshift slice to contain a high local num- 
ber density of galaxies (presumably many groups). There are 
more quiescent galaxies with low Sersic indices in this red- 
shift slice than are seen in other slices. This reinforces the 
notion that quiescent galaxies with low Sersic indices are the 
result of environment processes. 

If the environment of the cluster does suppress gas accre- 
tion onto galaxies, then we naturally expect a higher quies- 
cent fraction of galaxies in the vicinity of the cluster. We test 
this using the fraction of quiescent galaxies, /g, which is the 
ratio of the number of quiescent galaxies to the total number 
of galaxies. Taking galaxies with stellar mass > 2 x 10 10 M Q 
from the full UDS sample at the redshift of the cluster, we 
estimate an intrinsic quiescent fraction of /q = 0.53. This is 
consistent with the measurements from Quadri et al. (2012) 
considering differences in redshift binning, and in the mea- 
surements of the internal colors (see discussion in § 2, Quadri 
et al. 2012; Brammer et al. 2009). 

Figure 9 shows the fg as a function of projected distance 
from the center of the cluster. The figure also shows the prob- 
ability of obtaining the expected /q value for a binomial dis- 
tribution given the numbers of galaxies in each bin. Galaxies 
within 1 Mpc show an enhancement in the quiescent fraction, 
with /q=0. 62- 1.0. Using a binomial distribution, this is signif- 



icant with p = 0.033 for all galaxies within 1 Mpc, about the 
2cr level assuming a Gaussian distribution (and even higher 
significance, p = 0.01 1, considering the galaxies at 0.5-1 Mpc 
only). 

There is evidence for an enhanced quiescent fraction in the 
cluster. At projected distances /? pro j < 1 Mpc, we observe a 
significant enhancement (> 2ct) of the quiescent fraction (Fig- 
ure 9). In the annulus 1 Mpc < R pm j < 1.5 Mpc, the galaxies 
have a quiescent fraction that is consistent with the intrinsic 
fraction derived for all galaxies. As illustrated in Figure 4 the 
unchanging quiescent fraction from the annulus to the field 
results of a high number of star-forming galaxies in the annu- 
lus. We are unable to rule out the possibility that these star- 
forming galaxies result from contamination effects, and that 
the intrinsic quiescent fraction in the annulus is higher. It may 
also be possible that the some processes acting to suppress 
star-formation do not follow a spherically symmetric geome- 
try (which we assume using an annulus of constant projected 
radius). Larger samples of cluster galaxies observed at large 
projected distances are needed to make more rigorous conclu- 
sions. Nevertheless, the evidence here is that the galaxies in 
the 1 Mpc < R pro j < 1 .5 Mpc annulus have a quiescent fraction 
consistent with the field. This evidence disfavors (but does not 
exclude) the assembly bias scenario as this effect would likely 
exhibit an enhanced quiescent fraction in the near vicinity of 
the assembling cluster. 

Our findings at z = 1 .6 are similar to some studies at lower 
redshifts. Woo et al. (2012) show that at z ~ 0-0.2 and 
fixed stellar mass, the fraction of quiescent satellite galax- 
ies is strongly dependent on the projected distance from the 
center of the parent halo, and that the fraction of "quenched" 
galaxies is enhanced out to projected distances of ~ 1.5/? v ir 
around massive halos (M/, ~ 10 14 5 Mq). Weinmann et al. 
(2006) reach similar conclusions. Wetzel et al. (2011) simi- 
larly find that the fraction of quenched galaxies that are cen- 
trals in their halos is higher and deviates significantly from 
field values within 2 x R v i r , but at larger clustercentric dis- 
tances there is no indication that the environment has any ef- 
fect on the star-formation in galaxies. Depending on the size 
of the virial radius of the cluster (upper limit of /? v j r ~ 0.9 
Mpc, see above), it appears that the quiescent fraction is en- 
hanced out to distances of l/? v ir- Interestingly, the enhanced 
quiescent fraction does not extend to projected distances of 
1-1.5 Mpc, which would be expected if field galaxies are 
being quenched as they are accreted into the cluster environ- 
ment. A possible explanation is that the galaxies in the annu- 
lus are in a group sized dark matter halo(s), and may experi- 
ence excess star-formation (e.g., Tran et al. 2009) in addition 
to quenching that keeps the quiescent fraction approximately 
equivalent to the field value. 

Another interesting system is the well studied, X-ray lumi- 
nous cluster XMMU J2235. 3-2557 at a redshift of z = 1.39 
(Lidman et al. 2008; Strazzullo et al. 2010; Bauer et al. 
201 lb). The core of this highly evolved cluster exhibits a tight 
red sequence of galaxies (very similar to the cluster at z = 1 .62 
in our field) with little or no ongoing star formation. The clus- 
ter members on the outskirts on the other hand show greater 
diversity and a notable increase in [Oil] emitters and galaxies, 
which may be hosts of dust obscured star formation. There is 
also evidence for a "quenching radius" (located ~ 200 kpc 
from the cluster core) within which star formation is rapidly 
truncated, possibly in relation to the hot, X-ray emitting intr- 
acluster medium. Grutzbauch et al. (2012) extend the study 
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out to ~ 1 .5R V ir and find SFRs lower than typical values for 
field galaxies at the same redshift. It is noted that clusters de- 
tected due to their extended X-ray emission are more likely 
to be evolved structures than clusters detected by different 
methods. Observations of XMMU J2235. 3-2557 are consis- 
tent with more advanced stages of the evolution proposed for 
the cluster studied here. 

Therefore, our interpretation of the data is that environ- 
ments of higher density, such as the forming cluster at z = 
1.62, accelerate the morphological evolution and quiescent 
fraction of galaxies. Furthermore the data show that the en- 
vironment has two effects. First, Lotz et al. (2013) and Pa- 
povich et al. (2012) argue for an increased merger rate within 
the densest environments in order to account for the larger 
effective sizes of quiescent galaxies and the elevated number 
of faint companions observed. Second, here we argue that 
as galaxies enter the sphere of influence of the cluster, star- 
formation must be suppressed in some galaxies in order to 
explain their quiescent colors and morphologies (§ 4.4). The- 
ory predicts that the physical effects associated with the clus- 
ter could include gas-shock heating, as discussed above. An- 
other possibilty that can not be rejected is that galaxies near 
the cluster may experience assembly bias as they compete for 
baryonic gas with the larger gravitational potential of the clus- 
ter. Furthermore, Quadri et al. (2012) observe an upturn in 
the fraction of quiescent galaxies for lower-mass galaxies in 
higher density regions in the same mass and redshift range 
where we observe the change in morphologies, adding further 
evidence that the environment has a strong effect on galaxies 
of moderate mass. Therefore, at z = 1.6, it appears that the 
environment has a stronger effect on moderate-mass galaxies 
(2 x 10 10 -9 x 10 10 Mq), as it is these galaxies that exhibit the 
most difference in size and Sersic index compared to similar 
galaxies in the field (although we note that this applies less to 
the most massive quiescent galaxies, which show very similar 
sizes in both the cluster and field, Papovich et al. 2012). 

6. SUMMARY 

To summarize, in this paper we have studied the depen- 
dence of galaxy color, stellar mass, and morphology for galax- 
ies on environment at z = 1 .6 in the CANDELS/UDS field. 
This field containts Hubble Space Telescope near-infrared 
imaging over 9 '4 x 22 .'0, corresponding to a projected, phys- 
ical area of 4.8 Mpc x 11.2 Mpc at z = 1.6. This field also 
contains a known galaxy overdensity (cluster) at this redshift 
(Papovich et al. 2010). We define a sample of cluster galaxies 
as those within 1.5 Mpc of the cluster center, and we define a 
sample of galaxies in the field as those further than 3.0 Mpc 
of the cluster center. We also find it useful to define a sample 
of galaxies in an annulus between R pm j = 1-1.5 Mpc from the 
cluster center. Because we use the same suite of ground-based 
and space-based data to study the properties of the galaxies 
in our samples, our study is very homogenous. Any system- 
atic biases affecting our data or analysis will be identical, and 
studies of the relative properties between the cluster and field 
samples will be robust. 

We derive stellar masses, rest-frame colors, and SFRs using 
the broad suite of ground-based and Spitzer data. We quan- 
tify the morphology of the galaxies using GALFIT with the 
HST/F125W data to measure the effective radius, r e ff, and 
Sersic index, n, of all galaxies in our samples. In both the 
cluster and field, half of the bulge-dominated galaxies (« > 2) 
reside on the red sequence of the color-magnitude diagram, 
and most disk-dominated galaxies (n < 2) reside in the blue 



cloud associated with star-forming galaxies. Applying the 
WMW statistical test to field and cluster galaxies, we de- 
rive some evidence that the cluster galaxies have redder rest- 
frame U -B colors and higher stellar masses compared to the 
field, star-forming galaxies in both the cluster and field show 
no significant differences in their morphologies. In contrast, 
there is evidence that quiescent galaxies in the cluster have 
larger median effective radii, ^1" med = 2.0 kpc, and smaller 
Sersic indices, n^d" = 2.6 compared to the field, which have 
'■eff ld me d = l-4kpcand«S = 3.3. 

The differences in morphology are larger comparing quies- 
cent galaxies in an annulus defined by clustercentric distances 
1 Mpc < R pm j < 1.5 Mpc to the other quiescent galaxies. The 

quiescent galaxies in this annulus have median n l ^ d 5Mpc = 2.1 

an< ^ r reff med PC = 2.7 kpc, more consistent with the morpholo- 
gies of the star-forming galaxies. We find that the Sersic in- 
dex generally correlates with rest-frame U — V color, such that 
the majority of red, quiescent galaxies have higher Sersic in- 
dices (« > 2), consistent with the findings of Bell et al. (2012). 
However, a high fraction (w 50%) of the quiescent galaxies in 
the 1 Mpc < R pm j < 1 .5 Mpc annulus have n < 2 compared 
to only < 15% of quiescent galaxies overall. We argue that 
these galaxies have been processed under the influence of the 
cluster environment, and the evidence favors models where 
gas accretion onto these galaxies is suppressed. 

We argue that the primary channel for the evolution of 
galaxies at z = 1.6 depends on the mass of the galaxy, as the 
galaxies in both the cluster and field show weak differences in 
their masses, specific SFRs, and rest-frame U-B colors. 

Our results also show that the environment has some role 
in the morphological evolution of galaxies at z = 1.6. These 
differences are seen in the cluster as a lack of "compact" qui- 
escent galaxies observed in the field (see see also, Papovich 
et al. 2012). The quiescent cluster galaxies also have lower 
Sersic indices compared to the field, and this is especially true 
for the quiescent galaxies in the annulus of 1-1.5 Mpc from 
the cluster. The change in galaxy properties in this annulus 
is especially interesting as these clustercentric distances are 
greater than estimates of the virial radius, which have an upper 
limit of 0.9 Mpc from dynamical and X-ray constraints (see 
Papovich et al. 2010 and Pierre et al. 2011). Taken together, 
we argue that the environment influences star-forming galax- 
ies as they enter the sphere of influence of the cluster at ^2 
times the virial radius, and that this may affect more strongly 
the moderate mass galaxies. This is consistent with expecta- 
tions from theory, including effects from strangulation, pre- 
processing by groups, or assembly bias. Our data favors a 
scenario dominated by strangulation, however the other two 
possibilities cannot be excluded. 

As a final word, it is unknown if our results are extendable 
to all quiescent galaxies in different environments at z = 1 .6, 
including better statistics on group environments. Our results 
are based on a single forming galaxy cluster at z ~ 1 -6, and 
even with the large CANDELS dataset, we have morphologi- 
cal information from HST for fewer than 100 quiescent galax- 
ies at this redshift. Clearly larger samples of galaxies at this 
(and other) redshifts are required to generalize these results 
and to improve (or refute) the significance of the differences 
in morphology we observe here. Such tests will be possible 
using the full CANDELS dataset as well as other HST studies 
in fields of high redshift clusters. 



The Environmental Dependence of the Color-Mass-Morphology Relation at z = 1.6 



15 



We wish to thank the members of the CANDELS team for 
their contributions, and we thank L. Macri, and R. Quadri for 
stimulating discussions, and helpful comments. We also wish 
to thank the anonymous referee for comments and sugges- 
tions that improved the quality and clarity of this work. This 
work is based on observations taken by the CANDELS Multi- 
Cycle Treasury Program with the NASA/ESA HST, which is 
operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. This 
work is supported by HST program number GO- 12060. Sup- 
port for Program number GO- 12060 was provided by NASA 
through a grant from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Re- 
search in Astronomy, Incorporated, under NASA contract 
NAS5-26555. This work is based on observations made 
with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technol- 
ogy. This work is based in part on data obtained as part of 
the UKIRT Infrared Deep Sky Survey. MCC acknowledges 
support provided by NASA through Hubble Fellowship grant 
#HF-5 1269.01 -A, awarded by the Space Telescope Science 
Institute, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc., for NASA, under con- 
tract NAS 5-26555. MCC also acknowledges support from 
the Southern California Center for Galaxy Evolution, a multi- 
campus research program funded by the University of Califor- 
nia Office of Research. DC acknowledges receipt of a QEII 
Fellowship from the Australian Government. We acknowl- 
edge generous support from the George P. and Cynthia Woods 
Institute for Fundamental Physics and Astronomy. 

REFERENCES 



Abadi, M. G., Moore, B., & Bower, R. G. 1999, MNRAS, 308, 947 
Bahe, Y. M., McCarthy, I. G., Balogh, M. L., & Font, A. S. 2012, ArXiv 
e-prints 

Baldry, I. K., Balogh, M. L., Bower, R. G., et al. 2006, MNRAS, 373, 469 
Balogh, M. L., Navarro, J. F, & Morris, S. L. 2000, ApJ, 540, 1 13 
Barden, M., HauBler, B., Peng, C. Y, Mcintosh, D. H., & Guo, Y. 2012, 
MNRAS, 2576 

Bauer, A. E., Conselice, C. J., Perez-Gonzalez, P. G., et al. 2011a, MNRAS, 
417, 289 

Bauer, A. E., Griitzbauch, R., J0rgensen, I., Varela, J., & Bergmann, M. 

2011b, MNRAS, 411,2009 
Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 
Bell, E. F, van der Wei, A., Papovich, C, et al. 2012, ApJ, 753, 167 
Birnboim, Y, & Dekel, A. 2003, MNRAS, 345, 349 
Birnboim, Y, Dekel, A., & Neistein, E. 2007, MNRAS, 380, 339 
Blakeslee, J. P., Franx, M., Postman, M., et al. 2003, ApJ, 596, L143 
Blakeslee, J. P., Holden, B. P., Franx, M., et al. 2006, ApJ, 644, 30 
Blanton, M. R., Eisenstein, D., Hogg, D. W., Schlegel, D. J., & Brinkmann, 

J. 2005, ApJ, 629, 143 
Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503 
Brammer, G. B., Whitaker, K. E., van Dokkum, P. G., et al. 2009, ApJ, 706, 

L173 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Calzetti, D., Armus, L., Bohlin, R. C, et al. 2000, ApJ, 533, 682 
Chuter, R. W., Almaini, O., Hartley, W. G., et al. 2011, MNRAS, 413, 1678 
Cooper, M. C, Gallazzi, A., Newman, J. A., & Yan, R. 2009 
Cooper, M. C, Newman, J. A., Coil, A. L., et al. 2007, MNRAS, 376, 1445 
Cooper, M. C, Newman, J. A., Weiner, B. J., et al. 2008, MNRAS, 383, 
1058 

Cooper, M. C, Griffith, R. L., Newman, J. A., et al. 2011, MNRAS, 

submitted, arXiv: 1109.5698 
Croton, D. J., Gao, L., & White, S. D. M. 2007, MNRAS, 374, 1303 
Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 
Cucciati, O., Iovino, A., Kovac, K., et al. 2010, A&A, 524, A2 
Cuesta, A. J., Prada, F, Klypin, A., & Moles, M. 2008, MNRAS, 389, 385 
Daddi, E., Dickinson, M., Morrison, G., et al. 2007, ApJ, 670, 156 
Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 



Dekel, A., Birnboim, Y, Engel, G., et al. 2009, Nature, 457, 451 
Dressier, A. 1980, ApJ, 236, 351 

Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33 
Gao, L., Springel, V., & White, S. D. M. 2005, MNRAS, 363, L66 
Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35 
Griitzbauch, R., Bauer, A. E., J0rgensen, I., & Varela, J. 2012, MNRAS, 
423, 3652 

Griitzbauch, R., Conselice, C. J., Varela, J., et al. 2011a, MNRAS, 411, 929 
Griitzbauch, R., Conselice, C. J., Bauer, A. E., et al. 2011b, MNRAS, 418, 
938 

Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1 
Guo, Y, Mcintosh, D. H., Mo, H. J., et al. 2009, MNRAS, 398, 1 129 
Haas, M. R., Schaye, J., & Jeeson-Daniel, A. 201 1, MNRAS, in press, 
arXiv: 1103.0547 

Hansen, S. M., Sheldon, E. S., Wechsler, R. H., & Koester, B. P. 2009, ApJ, 
699, 1333 

Haussler, B., Mcintosh, D. H., Barden, M., et al. 2007, ApJS, 172, 615 
Hogg, D. W., Blanton, M. R., Eisenstein, D. J., et al. 2003, ApJ, 585, L5 
Hogg, D. W., Blanton, M. R., Brinchmann, J., et al. 2004, ApJ, 601, L29 
Holden, B. P., Blakeslee, J. P., Postman, M., et al. 2005, ApJ, 626, 809 
Homeier, N. L., Demarco, R., Rosati, P., et al. 2005, ApJ, 621, 651 
Kauffmann, G., White, S. D. M., Heckman, T. M., et al. 2004, MNRAS, 
353, 713 

Kocevski, D. D., Lemaux, B. C, Lubin, L. M., et al. 201 1, ApJ, 736, 38 
Koekemoer, A. M., Faber, S. M., Ferguson, H. C, et al. 2011, ApJS, 197, 36 
Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 
Kriek, M., van Dokkum, P. G., Franx, M., et al. 2006, ApJ, 649, L71 
Krist, J. 1995, in Astronomical Society of the Pacific Conference Series, 

Vol. 77, Astronomical Data Analysis Software and Systems IV, ed. 

R. A. Shaw, H. E. Payne, & J. J. E. Hayes, 349-+ 
Lidman, C, Rosati, P., Tanaka, M„ et al. 2008, A&A, 489, 981 
Lotz, J. M., Papovich, C, Faber, S. M., et al. 2013, ApJ, submitted, 

arXiv: 11 10.3821 

Mann, H. B., & Whitney, D. R. 1947, The Annals of Mathematical 
Statistics, 18, 50 

McGee, S. L., Balogh, M. L., Bower, R. G, Font, A. S., & McCarthy, I. G. 

2009, MNRAS, 400, 937 
Mei, S., Holden, B. P., Blakeslee, J. P., et al. 2009, ApJ, 690, 42 
Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A. 1996, Nature, 

379, 613 

Muldrew, S. I., Croton, D. J., Skibba, R. A., et al. 2012, MNRAS, 419, 2670 
Navarro, J. F. 1990, MNRAS, 242, 311 

Noeske, K. G., Weiner, B. J., Faber, S. M., et al. 2007, ApJ, 660, L43 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Papovich, C, Momcheva, I., Willmer, C. N. A., et al. 2010, ApJ, 716, 1503 

Papovich, C, Bassett, R., Lotz, J. M., et al. 2012, ApJ, 750, 93 

Patel, S. G., Holden, B. P., Kelson, D. D., Illingworth, G. D., & Franx, M. 

2009, ApJ, 705, L67 
Peng, C. Y, Ho, L. C, Impey, C. D., & Rix, H. 2002, AJ, 124, 266 
Peng, Y.-j., Lilly, S. J., Kovac, K, et al. 2010, ApJ, 721, 193 
Pierre, M., Clerc, N„ Maughan, B., et al. 2011, ArXiv e-prints 
Poggianti, B. M., Desai, V, Finn, R., et al. 2008, ApJ, 684, 888 
Postman, M., & Geller, M. J. 1984, ApJ, 281, 95 
Postman, M., Franx, M., Cross, N. J. G., et al. 2005, ApJ, 623, 721 
Quadri, R. F, Williams, R. J., Franx, M., & Hildebrandt, H. 2012, ApJ, 744, 

88 

Read, J. I., Wilkinson, M. I., Evans, N. W., Gilmore, G., & Kleyna, J. T. 

2006, MNRAS, 366, 429 
Rudnick, G. H., Tran, K.-V., Papovich, C, Momcheva, I., & Willmer, C. 

2012, ArXiv e-prints 
Santos, J. S., Rosati, P., Gobat, R., et al. 2009, A&A, 501, 49 
Sersic, J. L. 1968, Atlas de galaxias australes (Cordoba, Argentina: 

Observatorio Astronomico, 1968) 
Skibba, R. A., Bamford, S. P., Nichol, R. C, et al. 2009, MNRAS, 399, 966 
Strazzullo, V., Rosati, P., Pannella, M., et al. 2010, A&A, 524, A17 
Tasca, L. A. M., Kneib, J.-R, Iovino, A., et al. 2009, A&A, 503, 379 
Thomas, D., Maraston, C, Bender, R., & Mendes de Oliveira, C. 2005, ApJ, 

621, 673 

Tran, K, Saintonge, A., Moustakas, J., et al. 2009, ApJ, 705, 809 
Tran, K, van Dokkum, P., Franx, M., et al. 2005, ApJ, 627, L25 
Tran, K, Papovich, C, Saintonge, A., et al. 2010, ApJ, 719, L126 
van der Wei, A. 2008, ApJ, 675, L13 

van der Wei, A., Holden, B. P., Zirm, A. W., et al. 2008, ApJ, 688, 48 
van der Wei, A., Bell, E. F, Haussler, B., et al. 2012, ApJS, 203, 24 
Wechsler, R. H., Zentner, A. R., Bullock, J. S., Kravtsov, A. V., & Allgood, 
B. 2006, ApJ, 652,71 



16 



Bassett et al. 




21 22 23 24 25 26 

Output Magnitude 



Figure 10. The GALFIT output magnitude vs fractional difference in Sersic index (left panel) and effective radius (right panel) between input and output values. 
The blue points show the median values, the error bars depict the biweighted scale (Sbi) as described in Beers et al. (1990). The points are placed in half 
magnitude wide bins with the center of the first bin at 20.5. The panels are separated vertically by Sersic indices and effective radii values (respectively) used 
as input in our simulations. Note that because each panel is marginalized over the other parameter (e.g. n=l plot contains all r e /A this plot is not directly 
comparable to Table 2. Note also that the tendency for datapoints to be arranged in diagonal lines is due to the fact that input magnitudes were chosen at 0.5 
magnitude steps. This is similar to the simulations of Haussler et al. (2007). 
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APPENDIX 
GALFIT TESTS 

In order to determine to which magnitude we can reliably extract structural parameters of a galaxy using GALFIT, we per- 
formed a series of simulations. Galaxy models were created by running GALFIT with all fitted parameters (such as effective 
radius, Sersic index, and magnitude) fixed to chosen values. The models were then inserted randomly into the HST F125W 
image, avoiding the edge and extreme proximity or overlap with existing objects. Then GALFIT was run on these simulated 
galaxies embedded in our data in an attempt to recover the known values. The input models tested span a range of Sersic indices, 
effective radii, and magnitudes typical of galaxies in our sample. The output values from these tests were then separated into 
bins 0.5 mag in size and the robust median and biweighted scale (as described in Beers et al. 1990) were computed for various 
input parameters. The results are summarized in table 2, and plotted in Figure 10. As previously stated, our analysis of these 
simulations showed that the recovered effective radii are accurate to better than 40% and the Sersic indices to better than 25% for 
galaxies with r e g = 2 kpc and n = 4 with magnitude m(F125W) = 24 mag, near our stellar-mass limit of our sample, (2xlO lo M ). 
Galaxies of this type are among the most difficult to fit: the effective radius and Sersic index are highly covariant. Uncertainties 
are significantly lower for more disk-like and brighter sources (as shown in Table 2). 
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Table 2 

Summary of GALFIT tests for HST WFC3 F125W Data 



Output Magnitude r out Range (kpc) 


n om Range 


■ 

SBi{(r ol u-r in )/r in ] 




SBi[{noui-ni„)/n in \ 


2().5±0.25 0-.75 


<1 


0.002 


0.005 




1-2.5 


0.003 


0.010 




>2.5 


0.010 


0.016 


1.0±.25 


<1 


0.002 


0.004 




1-2.5 


0.006 


0.009 




>2.5 


0.018 


0.018 


1.5±.25 


<1 


0.004 


0.005 




1-2.5 


0.009 


0.012 




>2.5 


0.024 


0.020 


2.0±.25 


<1 


0.005 


0.005 




1-2.5 


0.012 


0.013 




>2.5 


0.031 


0.021 


21.0±0.25 0-.75 


<1 


0.003 


0.008 




1-2.5 


0.005 


0.018 




>2.5 


0.017 


0.026 


1.0±.25 


<1 


0.004 


0.007 




1-2.5 


0.009 


0.015 




>2.5 


0.031 


0.031 


1.5±.25 


<1 


0.005 


0.008 




1-2.5 


0.015 


0.020 




>2.5 


0.044 


0.034 


2.0±.25 


<1 


0.008 


0.010 




1-2.5 


0.019 


0.020 




>2.5 


0.065 


0.044 


21.5±0.25 0-.75 


<1 


0.005 


0.012 




1-2.5 


0.009 


0.026 




>2.5 


0.025 


0.041 


1.0±.25 


<1 


0.006 


0.010 




1-2.5 


0.017 


0.028 




>2.5 


0.046 


0.046 


1.5±.25 


<1 


0.010 


0.015 




1-2.5 


0.025 


0.033 




>2.5 


0.067 


0.049 


2.0±.25 


<1 


0.013 


0.015 




1-2.5 


0.031 


0.035 




>2.5 


0.101 


0.068 


99 n+C) 25 0- 75 


<J 


0.006 


0.017 




1-2.5 


0.015 


0.042 




>2.5 


0.048 


0.078 


1.0±.25 


<1 


0.011 


0.017 




1-2.5 


0.025 


0.040 




>2.5 


0.059 


0.065 


1.5±.25 


<1 


0.015 


0.023 




1-2.5 


0.036 


0.048 




>2.5 


0.086 


0.071 


2.0±.25 


<1 


0.017 


0.027 




1-2.5 


0.053 


0.058 




>2.5 


0.256 


0.158 


22.5±0.25 0-.75 


<1 


0.014 


0.030 




1-2.5 


0.024 


0.065 




>2.5 


0.068 


0.109 


1.0±.25 


<1 


0.018 


0.022 




1-2.5 


0.042 


0.065 




>2.5 


0.113 


0.106 


1.5±.25 


<1 


0.023 


0.034 




1-2.5 


0.056 


0.070 




>2.5 


0.140 


0.120 


2.0±.25 


<1 


0.031 


0.047 




1-2.5 


0.087 


0.090 




>2.5 


1.033 


0.643 


23.0±0.25 0-.75 


<1 


0.018 


0.047 




1-2.5 


0.038 


0.092 




>2.5 


0.095 


0.161 


1.0±.25 


<1 


0.027 


0.039 




1-2.5 


0.063 


0.106 




>2.5 


0.152 


0.161 


1.5±.25 


<1 


0.036 


0.051 




1-2.5 


0.089 


0.110 




>2.5 


0.222 


0.196 


2.0±.25 


<1 


0.053 


0.070 




1-2.5 


0.129 


0.135 




>2.5 


3.571 


0.881 


23.5±0.25 0-.75 


<1 


0.033 


0.074 




1-2.5 


0.056 


0.154 
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Table 2 — Continued 



Output Magnitude r out Range (kpc) n out Range S B i[(r„ u , -r in )/r in ] S B i[{n„„, -n in )/n in \ 





>2.5 


0.153 


0.253 


1.0±.25 


<1 


0.044 


0.066 




1-2.5 


0.089 


0.148 




>2.5 


0.213 


0.223 


1.5±.25 


<1 


0.063 


0.074 




1-2.5 


0.116 


0.161 




>2.5 


0.275 


0.231 


2.0±.25 


<1 


0.093 


0.094 




1-2.5 


0.200 


0.202 




>2.5 


3.360 


1.019 


24.0±0.25 0-.75 


<1 


0.060 


0.132 




1-2.5 


0.108 


0.273 




>2.5 


0.186 


0.301 


1.0±.25 


<1 


0.071 


0.098 




1-2.5 


0.138 


0.251 




>2.5 


0.332 


0.308 


1.5±.25 


<1 


0.086 


0.123 




1-2.5 


0.194 


0.285 




>2.5 


0.384 


0.359 


2.0±.25 


<1 


0.132 


0.189 




1-2.5 


0.273 


0.305 




>2.5 


3.243 


0.940 
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